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ABSTRACT 


Arrays  of  performance  value*  shoving  the  change  in  noise 
variance  due  to  degho.ting  indicate  significant  sensitivity  to 
both  the  deghost ing  method  and  pre-band  pass  filtering.  As  a 
specification  for  single  channel  deghosting,  the  apparent  re¬ 
flection  coefficient  used  to  derive  the  operator,  should  be 
constrained  such  that  the  noise  amplification  does  not  exceed 
specified  design  criteria  such  as  1  db.  Some  typical  results 
of  such  a  test  are  ?iven.  The  apparent  reflection  coefficient 
used  to  design  the  echo  suppressor  is  observed  as  highly  depen 
dent  on  the  echo  time. 


1 •  INTRODUCTION 

Single  channel  deghoating  is  considered  in  the  context  of 
a  modular  approach  to  processing  vertical  array  data  for  signal 
detection  and  measurement.  The  modules  consist  oft 

(1)  pre-set  band  pass  filter  designed  to  improve  detection  pro¬ 
bability  of  nearly  vertical  p-wavee  for  a  fixed  false  alarm  pro¬ 
bability  (eg  0.02  or  expectation  of  less  then  one  false  alarm  per 
minute ) 

(2)  single  channel  deghosting  operator  for  removal  of  echo  dis¬ 
tortion 

(3)  possible  multi-channel  wave  number- frequency  filter,  fan 
filter,  for  estimating  the  p-pulse  and  suppressing  the  ambient  and 
signal-coda  noise 

(4)  in  place  of  (3),  possible  use  of  signal  correlation  criteria 
for  detecting  and  measuring  the  p-pulse. 

Here,  we  are  confined  to  evaluating  the  performance  of  single 
channel  degnosting  operators  to  determine  the  constraints  required 
to  prevent  the  deghosting  operatorr  from  degrading  the  overall 
performance  of  the  system. 

I .  NOISE  DATA 

The  noise  j.s  obtained  from  a  pseudo  random  number  generating 
function,  FORTRAN  63  RANF  (-1).  The  alogorithm  has  been  tested  and 
shown  statistically  to  provide  a  sequence  of  random  uncorrelated 
numbers.  The  sequence  will  not  repeat  until  millions  of  samples 
have  been  called.  The  pseudo  random  numbers  are  averaged  nine  times 
to  produce  a  nearly  Gaussian  input  of  uncorrelated  samples.  They 
are  filtered  twice  in  series  using  RECFIL3,  a  recursive  narrow-band 
filter.  The  low  frequency  noise  component  is  obtained  with  the 
frequency  set  at  to  .2  cps  and  Q  =  fo/Af  *  3.5,  and  the  high  fre¬ 
quency  to  2.0  cps,  Q  ■  6.  A  plot  of  the  power  spectral  density 
of  a  sample  of  the  synthetic  noise  is  shown  in  figure  1. 


3 •  PERFORMANCE  ANALYSIS  PARAMETERS 

The  purpose  of  applying  a  single  channel  inverse  deghosting 
operator  is  to  suppress  the  echo  reflected  at  the  earth's  surface. 

For  clcsfisinit  signals,  the  Incidence  angle  is  approximately  15°, 
thus  the  down-going  pulse  could  bo  modeled  with  a  reflection  coef- 
ficent  of  0.9  and  the  up-going  pulse  with  a  reflection  coefficient 
of  1.1.  in  practice  however,  reflection  coefficients  greater  than 
1.0  leal  to  unstable  operators,  so  that  single  channel  echo-suppress¬ 
ion  is  feasible  only  if  we  attempt  to  only  partially  remove  the  echo 
in  isolating  the  up-going  and  down-going  pulses.  In  analysing  the 
performance  of  our  echo-suppression  filters,  we  try  different  reflect¬ 
ion  coefficients  from  zero  to  .9  at  intervale  of  approximately  .1, 
further  w*  sample  two-way  echo  times  from  .05  to  1.2  seconds  at 
intervals  ->f  0.5  seconds.  Four  minute  noise  samples  are  filtered. 

For  each  sample  specified  by  the  arrays  of  echo  times  and  reflection 
coefficients,  the  variance  of  the  noise  is  computed  before  and  after 
applying  the  deghosting  operator.  The  performance  of  the  filter  is 
guaged  in  decibels  as  10  log  (variance  after  filtering Aarianc. 
before  filtering).  Negative  values  indicate  a  reduction  in  noise; 

positive  an  amplification  in  noise  due  to  applying  the  degho.ting 
filter. 

^ '  pre_-filtering  and  deghosting  models 

This  analysis  will  demonstrate  that  the  performance  on  ambient 
noise  is  very  sensitive  to  the  bandpass  filter  and  also  on  the  num¬ 
erical  method  used  to  deghost.  Two  different  deghosting  methods,  two 
different  bandpass  filters,  and  unfiltered  data  make  up  a  total  of 
five  cases  to  be  evaluated  by  performance  arrays.  To  produce  each 
element  in  the  array  we  operate  on  four  minutes  of  data:  each  array, 

16  hours  of  data,  and  for  all  five  cases,  80  hours  of  data.  The 
fact  that  all  of  these  performances  figures  were  computed  with  less 
than  40  minutes  of  computer  time  attests  to  speed  of  the  bandpass 
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and  degho.ting  operator*. 

The  degho.ting  method.  are  described  am  method  I  and  method  II. 
Both  are  recursive  filter,  which  require  only  one  multiplication 
per  data  point . 


Method  I 


Method  II 


X  x  -  R 

r  ♦  m'2 


r-m/2 


mi  ♦  J 


ra(  i- 1  )♦ j  Ymi  +  j 
M  <  Ml  <  n 
1  <  j  <  M-l 


xr  i.  the  array  of  data  containing  the  echo 
R,  reflection  coefficient 
P  »  degho.tcd  data 

m  >  echo-time  measured  in  number  of  point. 


mi  ♦  j 


r  +  m/2 
R  •  p 


0.9 


r-m/2 


»(i-l)+j  +  Yni  +  j 


The  results  evaluated  for  two  band-pa.,  filter,  employ  gapped 

finite  difference  which  involve,  only  several  subtraction,  per  data 
point . 

5.  GAPPED  FINITE  DIFFERENCE  PTI/rER 

■  f  y 

The  gapped  difference  operator  is  given  a.  \  1,  o.o.  .  .o.  -k/, 

where  k  <  1,  and  the  qap  i.  specified  by  the  number  of  point,  between 

1  and  k.  Here  we  are  given  by  k  -  1 .  The  spectrum  of  the  operator 
is  given  by 


D(f ) 


i 


f  "s 

fo  J 


where  the  frequency  of  the  first  peak  of  the  filter  is  related  to 

the  number  of  points  in  the  gap,  L. 

S.R 

~o 


-  3  - 


L  * 


8.  R.,  the  sampling  rate. 


Note  that  series  application  of  D(f)  can  be  used  to  generate  higher 
differences  m 

D(m)  (fo)  -  a2)  7  -  <-nN  -d  n  <fQ> 

where  n  -  1,2,3.  ...  specifies  the  2nd  difference,  4th  difference, 

etc. Taken  as  a  filter  we  observe  that  for  n  *  odd,  the  sign  of  the 
filter  is  reversed.  The  response  of  the  filter  is 

(f)  -  (-I)"  «in2n  1  2  fo  J 

The  gap,  L,  is  determined  by  sett  mg  the  firct  peak  to  the  desired 
frequency.  Nulls  occ**r  at  integer  multiples  of  twice  the  frequency 
of  the  first  peak,  and  additional  peaks  occur  at  frequencies  between 
the  null's.  Suppose  at  frequencies  greater  than  the  first  null  in 
the  spectrum,  for  example  at  f-2,  the  aigral  and  noise  power  are 
expected  to  attenuate  rapidly  with  increasing  frequency.  These 
high  frequency  peaks  will  in  this  case  produce  negligible  effects 
in  the  filtered  output.  If,  however,  we  wish  to  insure  attenuation 
of  high  frequency  peaks  in  the  noise  spectrum  wg  may  apply  in  series 
with  the  filter  a  low  pass  sracothirg  operator  ^1.  o,  o,  o,  +kj 

where  k  <  1.  The  spectrum  of  the  operator  is  given  by  S(f). 

S(f)  -  cos  (-  y  -~) 

Repeated  application  of  the  smoothing  operator  leads  to  the  operator 


where,  for  example,  fo  may  be  taken  at  the  toldirg  frequency. 

For  k  -  1  there  are  no  mult iplicat lone  in  either  the  smoothing 
or  difference  operator,  so  chat  the  above  digital  operators  are 
exceedingly  fast.  For  k  -  1  only  an  approximate  inverse  can  be  ob¬ 
tained.  These  *re,  in  f'iCt,  simply,  deghosting  operators.  k  <  1, 

exact  inverses  can  be  retained .  The  invars*  can  b"  designed  as  re¬ 
cursive  filters  with  n  multiplications  per  point  where  n  is  the  number 


(n)  c( n) 

»f  series  applications  of  r  or  » 

Smoothed  gapped  finite  differenoe  filter,  can  he  de.cribed 


spectrally 


F  ( f )  -  D 


(2n)  (f)  s(,n)  (f) 


n  2n  1  tt  f  \ 

(-D  Sin  V - fT"  J 


The  two  models  evaluated  were 


1)  Smoothed  second  difference  filter 

n  -  1  fQ  -  1.1. cps  m  -  4  fx  -  10  cps 

2)  Smoothed  fourth  difference  filter 

n  =  2  f  ■  1  1  cps  m  =  4  fl  "  10‘cps 

6.  RESULTS 

Figure  one  .how.  the  power  spectrum  of  the  noise  data  used  in 
the  evaluation  of  combined  bandpass  filters  and  deghosting  operator.. 
Table  I  through  Table  V,  for  different  values  of  the  reflection  coe¬ 
fficient  and  the  echo  time  (sum  of  up-hole  and  down-hole  time)..  ow 
the  effect  of  the  deghosting  on  the  noise.  Table  1  is  for  un 
data,  on  the  other  tables  the  noise  is  bandpass  filtered.  The  per¬ 
formance.  figure,  in  decibel,  .how  the  effect  of  degho.ting  by  com- 
paring  variance  of  the  filtered  noise  after  deghosting  to  the  variance 
o£  filtered  noise  before  deghosting.  Examination  o,  the  performance 
figures  show  that  both  the  bandpass  filter  and  the  deghosting  method 
significantly  affect  the  performance  of  single  channel  deghosting 

operators . 

7.  CPU  -UlSIOHS 

For  unfiltered  data  and  echo  times  between  .2  sec  and  0.35 
seconds,  and  greater  than  0.6  seconds  the  assumed  reflection  coeffi¬ 
cient  should  not  exceed  0.75  which  is  adequate  for  removal  of 
of  the  echo.  For  data  prefiltered  by  smoothed  second 

operators,  the  best  performance  is  obtained  with  degho.ting  method  I. 
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For  echo  times  greater  than  0.3  seconds,  the  reflection  coefficient 
should  not  exceed  0.4  for  removal  of  slightly  less  than  half  of  the 
echo  amplitude. 

For  data  pretiltered  by  fourth  difference  operators,  the  best 
performance  was  obtained  usinc  deghosting  method  I.  For  echo  times 
exceeding  0.3  seconds  and  less  than  0.75  seconds,  the  reflection 
coefficient  should  rot  exceed  0.2  seconds.  For  echo  tiroes  greater 
than  0.75  but  less  than  0.9  seconds,  the  reflection  coefficient 
should  not  exceed  0.45;  for  echo  times  equal  or  greater  than  0.9 
seconds  the  reflection  coefficient  should  not  exceed  0.75.  The 
fourth  difference  filter  appears  most  suited  as  the  bandpass  filter. 
Due  to  deghosting,  there  will  be  little  if  any  compromise  of  the 
expected  9  db  gain  obtained  by  applying  the  bandpass  filter.  Any 
residual  echo  distortion  must  be  removed  by  multichannel  filtering. 

For  display  purposes,  an  adequate  job  of  deghosting  can  be 
done  on  all  channels  in  the  absence  of  bandpass  filtering.  Provided 
that  the  assumed  refl  ■rtion  coefficient  does  not  exceed  .75,  less 
than  1  db  total  noise  amplification  is  assured. 


RELATIVE  POWER 


Figure  1.  Noise  Power  Spectrum 
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Smoothed  Second  Difference 


Table  3.  Smoothed  Fourth  Difference  Filter,  Dtghosting  Method  I. 
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Table  5.  Smoothed  Fourth  Difference  Filter,  Deghosting  Method  II. 
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13  ABSTRACT 

Arrays  of  performance  values  showing  the  change  in  additive 
noise  variance  due  to  single  channel  deghosting  indicate  signifi¬ 
cant  sensitivity  to  both  the  deghosting  method  and  band  pass 
filtering.  Only  part  of  the  echo  can  be  removed  feasibly,  for 
example,  with  less  than  one  db  increase  in  the  noise,  by  the 
single  channel  deghosting  filter.  The  maximum  fraction  of  the 
echo  which  can  be  removed  varies  greatly  with  the  echo  time, 
especially  if  the  noise  was  bana  pass  filtered. 
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